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ABSTRACT 


The  shear  in  the  mixed  layer  was  forecasted  and  studied  for  the  duration  of 
the  MILE  experiment.  The  forecasting  was  done  with  the  Warn-Varnas  and  Piacsek 
(1979)  model. 

The  forecasts  were  initiated  from  data  and  driven  by  the  experimentally- 
measured  wind  stress,  radiation,  and  latent  and  sensible  heat  flux.  The  duration 
of  forecasts  spanned  a  time  frame  of  several  days  to  a  month. 

The  problem  of  understanding  the  shear  in  the  mixed-layer  region  during  MILE 
was  approached  from  two  points  of  view.  First,  the  basic  dynamic  processes  that 
generate  shear  were  considered.  Second,  statistical  distributions  of  shear  as  a 
function  of  depth  and  time  were  computed. 

The  basic  dynamics  can  be  understood  by  decomposing  the  velocity  into  an 
Ekman-type  component  and  an  inertial  component.  The  vertical  structure  of  this 
velocity  can  be  visualized  as  a  superposition  of  an  inertial  oscillation  onto  an 
Ekman  spiral;  the  inertial  oscillations  being  excited  by  abrupt  changes  in  wind 
stress.  In  most  cases,  the  maximum  shear  occurred  at  the  bottom  of  the  mixed 
layer  and  near  the  surface.  In  some  cases,  two  shear  maxima  regions  were  found 
below  the  surface.  These  cases  involved  situations  of  rising  wind  stress  from 
a  previous  history  of  low  magnitude. 

In  the  statistical  approach,  scatter  diagrams  of  shear  versus  Brunt-Vaisala 
frequency  are  computed.  At  different  depths,  these  scatter  diagrams  show  the  re¬ 
lationship  of  mixing,  as  indicated  by  the  Brunt-Vaisala  frequency,  to  shear.  The 
simulated  magnitudes  of  shear  are  in  the  range  of  the  experimental  results  of 
R.  E.  Davis ,  et  al . 

Probability  diagrams  of  calculated  flux  Richardson  numbers  as  a  function  of 
depth  are  also  computed.  These  probability  diagrams  show  a  shift  toward  higher 
Richardson  numbers  as  the  depth  increases.  This  reflects  the  fact  that  the  pene¬ 
tration  of  turbulence  and  mixing  tend  to  decrease,  on  the  average,  with  depth. 
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AN  ANALYSIS  OF  MODELLED  SHEAR  DISTRIBUTIONS  DURING 
THE  MILE  EXPERIMENT 

I.  INTRODUCTION 

The  propagation  of  internal  waves  in  the  ocean  is  strongly  affected  by  gradi¬ 
ents  in  density  and  current  velocity.  Such  regions  cause  reflection,  refraction, 
absorption,  and  breakinq  of  the  waves.  Under  certain  conditions,  the  waves  can 
also  grow  in  regions  of  velocity  shear.  To  be  able  to  predict  the  propagation  of 
a  group  of  internal  waves  in  the  upper  ocean,  one  must  be  able  to  predict  in  some 
deterministic/statistical  sense  the  properties  of  this  upper  ocean  environment;  in 
particular,  the  spatial  and  temporal  distribution  of  density  and  velocity  graidents 
and  regions  of  turbulence  activity.  Since  temperature  measurements  are  much  easier 
to  carry  out  than  current  measurements ,  and  temperature  data  sets  are  in  great 
abundance,  emphasis  will  be  on  measurements  and  theoretical  understanding  of  cur¬ 
rent  shear  in  the  coming  year. 

The  present  study  deals  with  the  problem  of  current  shear  in  the  mixed  layer. 
The  shear  in  the  mixed  layer  was  forecasted  and  studied  for  the  duration  of  the 
MILE  experiment.  The  forecasting  was  done  with  the  Warn-Varnas  and  Piacsek  (1979) 
model ,  referred  hereafter  as  the  WP  model . 

The  WP  model  is  a  one-dimensional,  higher-order  closure  model  in  the  sense 
that  it  retains  the  tri ple-correl ation  terms.  The  model  contains  forecast  equa¬ 
tions  for  the  mean  fields,  second-order  correlations,  and  third-order  correlations. 
The  finite-difference  solution  of  the  equations  is  formulated  without  any  critical 
Richardson  number  cut-off  criteria  of  the  turbulence,  and,  therefore,  the  effects 
of  shear  penetration  of  turbulence  at  the  bottom  of  the  mixed  layer  are  preserved. 

The  forecasts  were  initiated  from  data  and  driven  by  the  experimentally- 
measured  wind  stress,  radiation,  and  latent  and  sensible  heat  flux.  The  duration 
of  forecasts  ranged  from  several  days  to  a  month. 

The  problem  of  understanding  the  shear  in  the  mixed  layer  during  MILE  was 
approached  from  two  points  of  view.  First,  the  basic  dynamic  processes  that  gen¬ 
erate  shear  were  considered.  Second,  the  statistical  distributions  of  shear  as  a 
function  of  depth  and  time  were  calculated. 

The  basic  dynamics  can  be  understood  by  decomposing  the  velocity  into  an 
Ekman-type  component  and  an  inertial  component.  The  Ekman  components  are  obtained 
from  a  balance  between  the  second-order  correlation  terms  and  the  Coriolis  terms. 
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The  vertical  structure  of  this  velocity  can  be  visualized  as  a  superposition  of  an 
inertial  oscillation  onto  an  Ekman  spiral,  the  inertial  oscillations  being  excited 
by  abrupt  changes  in  wind  stress. 

This  hypothesis  was  illustrated  by  turning  on  a  constant  wind  stress,  with 
no  heat  flux,  for  one  of  the  MILE  temperature  profiles  and  performing  a  numerical 
simulation  of  several  inertial  periods.  Throughout  the  mixed  layer,  the  velocity 
hodographs  showed  circles  of  constant  amplitudes  versus  depth.  The  center  of  the 
circles  versus  depth  were  displaced  as  if  sitting  on  top  of  an  Ekman  spiral.  The 
amplitude  of  this  inertial  velocity  was  constant  throughout  the  mixed-layer  depth 
and  was  zero  below  the  mixed  layer.  The  maximum  shears  occurred  at  the  bottom  of 
the  mixed  layer  and  near  the  surface.  Both  components  of  the  velocity,  the  Ekman 
and  inertial,  showed  appreciable  shear  at  the  bottom  of  the  mixed  layer. 

Next,  a  time  was  selected  during  MILE  when  there  was  no  wind  stress.  In 

this  case,  it  was  found  that  the  Ekman  part  of  the  velocity  was  zero  and  the  shear 

at  the  bottom  of  the  mixed  layer  was  due  to  the  inertial  component  of  the  velocity 

only.  In  another  case,  a  time  was  selected  when  the  stress  peaked  around  4.5 
2  -2 

dynes/cm  ,  and  a  shear  of  5  x  10  cm/sec  was  found  at  the  bottom  of  the  mixed  lay¬ 
er.  The  magnitudes  of  the  Ekman  versus  the  inertial  parts  of  the  shear  were  3  x 
-  2  -6 

10  cm/sec  and  4  x  10  cm/sec,  respectively.  Also  considered  was  a  case  that  in¬ 
volved  a  situation  where  the  wind  stress  was  rising  from  a  previous  history  of  low 
magnitude.  The  maximum  shear  occurred  at  the  penetration  depths  of  the  turbulence 
generated  by  this  rising  wind  stress.  The  next  largest  shear  was  located  further 
down  at  the  bottom  of  the  mixed  layer  and  was  due  to  a  previously-excited  inertial 
wave. 

In  the  statistical  approach,  scatter  diaqrams  of  shear  versus  Brunt-Vaisal a 
frequency  are  computed.  At  different  depths,  these  scatter  diagrams  show  the 
relationship  of  mixing  to  shear  as  indicated  by  the  Brunt-Vai sa 1  a  frequency.  The 
simulated  magnitudes  of  shear  are  in  the  range  of  the  experimental  results  of 
R.  E.  Davis ,  et  al . 

Probability  diagrams  of  calculated  Richardson  numbers  as  a  function  of  depth 
are  also  computed.  These  probability  diagrams  show  a  shift  toward  higher  Richard¬ 
son  numbers  as  the  depth  increases.  This  reflects  the  fact  that,  the  penetration 
of  turbulence  and  mixing  tend  to  decrease,  on  the  average,  with  depth. 
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II.  EQUATIONS  OF  THE  ONE-DIMENSIONAL  WP  MODEL 
A.  The  Mean  Field  Equations 
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where  pQ  is  the  mean  density,  F  is  the  downward  flux  of  solar  radiation, 
is  the  specific  heat  at  constant  pressure. 

B.  The  WP  Equations  for  the  Second-Order  Correlations 
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where  e  is  the  dissipation  function. 
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The  constants  are  set  according  to  Mellor  (1973);  =  3.19,  C2  =  0.056, 

C.,  -  1.47,  Clt  =  4.78,  and  C4  =  1.87.  They  are  based  mostly  on  experimental  mea¬ 
surements  in  the  neutral ,  constant -flux  layer  and  grid-generated,  constant  shear 
+  l  nw . 

The  equations  for  the  triple  correlations  are  given  in  the  WP  paper  with  the 
decay  terms  changed  to  the  form  -2.48  x  triple  correlations.  The  finite- 
difference  scheme  used  to  solve  the  equations  is  described  in  the  Appendix. 

C.  Boundary  Conditions 

At  the  bottom  boundary,  all  velocities  are  set  to  zero  and  the  temper¬ 
ature  is  prescribed.  At  the  top  boundary,  or  surface,  one  uses 

“-V  =  ■  \ 


wT  -  -  H 


where  H  is  the  sum  of  the  latent,  sensible,  and  back  radiation  fluxes,  and  and 
r  are  the  wind  stresses  along  the  x  and  y  axis. 

D.  Dissipation  Function 

The  dissipation  function  is  defined  as 

r  k  1 ' /\  (2.8) 

-  -  .1886  as  obtainr-'1  from  Mellor  (  1973).  For  the  WP  model,  ■  is  defined 


(2.9) 
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is  Von  Karman’s  constant,  and  m,  is  an  empirical  con- 
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III.  RESULTS  AND  DISCUSSION 

A.  Statistical  Distributions  of  Stratifications  and  Shear 

A  numerical  simulation  of  the  MILE  experiment  with  the  WP  one-dimen¬ 
sional  model  was  performed.  The  data  used  was  that  of  Miyake  (1978).  The  driving 
fluxes  that  arose  from  this  data  are  shown  in  figure  1.  The  model  was  initialized 
from  measured  temperature  profiles  at  10:00z  August  2,  1977.  Then,  the  measured 
wind  stress,  solar  radiation,  back  radiation,  and  sensible  and  latent  heat  flux 
forced  the  model  via  the  boundary  conditions  as  given  in  Warn-Varnas  et  al . 

(1981).  Probability  distributions  of  flux  Richardson  numbers  were  calculated  at 
various  depths  for  the  duration  of  MILE  (about  30  days).  The  results  are  shown 
in  figures  2  through  4.  The  flux  Richardson  number  was  defined  as 


-ag  w'T' 

RF  =  - r - -  ’  (3-D 

.  7V  — 

dZ  *Z 

and  the  observations  were  4  hours  apart.  Near  the  surface,  inversions  in  temper¬ 
ature  can  occur  due  to  cooling.  This  results  in  negative  values  of  the  flux 
Richardson  number.  In  general,  the  probability  of  finding  the  low-value  Richard¬ 
son  numbers  is  greatest  near  the  surface  since  the  chances  of  turbulence  existing 
there  are  the  greatest.  For  turbulence  to  penetrate  further  into  the  mixed  layer, 
larger  wind  stresses  or  cooling  fluxes  are  required.  These  larger  values  of  forc¬ 
ing  fluxes  are  less  probable,  and,  therefore,  the  values  of  the  Richardson  numbers 
increase  with  depth,  indicating  that  there  is  less  turbulence  at  greater  depths. 
This  fact  is  illustrated  by  the  shift  of  the  probability  distribution  diagrams 
toward  higher  Richardson  numbers  as  the  depth  increases  (figures  2  through  4). 

An  estimate  of  the  distribution  of  stratification  and  shear  during 
MILE  can  be  obtained  from  scatter  diagrams  for  various  depths.  These  diagrams 
have  been  computed  from  the  model  results  at  intervals  of  4  hours.  The  results 
are  shown  in  figures  5  through  7.  The  units  are  cgs. 

In  general,  the  largest  shears  are  around  the  bottom  of  the  mixed 
layer,  which  is  located  around  30  meters  most  of  the  time  (figure  7).  The  next 
largest  shears  are  located  around  the  surface  (figure  5).  The  region  which  is 
the  mixed  layer  most  of  the  time  shows  smaller  shears  since  the  velocities  there 
change  less  with  depth  (figure  6).  The  largest  strati fications  are  near  the 
bottom  of  the  mixed  layer  where  the  thermocline  is  located  (figure  7). 
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B.  Case  of  Constant  Wind  Stress  With  No  Heat  Flux 


In  order  to  illustrate  the  consequences  of  wind-generated  shear  turbu- 

2 

lence  on  mixed-layer  dynamics,  a  constant  wind  stress  of  Tx  =  Ty  =  dynes/cm 
is  prescribed,  and  the  problem  is  initialized  from  the  MILE  temperature  profile  on 
day  number  7.  The  results,  about  four  inertial  periods  later,  are  shown  in  fig¬ 


ures  8  through  13.  The  velocity  can  be  decomposed  into  Ekman  and  inertial  compo¬ 
nents  as  v  =  v^  +  Vj.  It  is  assumed  that  the  Ekman  or  shear-type  components  of 
the  velocity  are  determined  by  the  equations 


u 


1  a  v 1  w 1 
E  "  f  az 


(3.2) 


Nomitsu  (1933)  gives  further  details  on  an  analytical  analog  of  such 
a  decomposition. 

The  largest  shear  is  located  near  the  bottom  of  the  mixed  layer  (fig¬ 
ure  11).  There,  the  Ekman  and  inertial  components  of  the  velocity  both  contribute 
to  the  shear  (figures  9  and  10).  An  interesting  result  is  shown  in  figure  10;  the 
amplitude  of  the  inertial  velocity  is  constant  with  depth  and  the  mixed  layer 
oscillates  as  a  "slab" . 


Figures  12  and  13  show  velocity  hodographs  at  depths  of  1  and  9  meters. 
Again,  one  can  see  a  constant  amplitude  oscillation  in  the  mixed  layer.  A  quali¬ 
tative  picture  of  the  dynamics  can  be  constructed  by  visualizing  an  Ekman  spiral 
oriented  about  50  degrees  to  the  right  of  the  wind  stress  with  a  constant  ampli¬ 
tude  inertial  oscillation  superposed  on  it.  Thus,  as  one  goes  down  into  the  ocean, 
the  hodographs  shift  to  the  right  just  as  if  they  were  following  the  Ekman  spiral. 

C.  Low  and  High  Wind  Stress  Cases  During  MILE 

Around  day  number  8,  the  wind  stress  peaks  at  about  4  dynes/cm  and, 
then,  around  day  8.5  diminishes  to  zero  (figure  1).  The  results  for  day  number  8 
are  shown  in  figures  14  through  16.  As  in  the  previous  example  of  constant  wind 
stresses,  the  Ekman  and  inertial  parts  of  the  velocity  solution  contribute  to  the 
shear  at  the  bottom  of  the  mixed  layer.  The  resultant  shear  magnitude  is  about 
0.1  cm/sec  (figure  16).  On  day  number  8.58,  both  the  wind  stress  and  the  Ekman 
part  of  the  velocity  solution  tend  toward  zero.  The  total  velocity  is  now  equal 
to  the  inertial  component  of  the  velocity  shown  in  figure  17.  The  resultant  'hear 
is  shown  in  figure  18.  Thus,  even  in  a  case  when  there  is  no  wind  stress  and  no 
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Ekman  component  of  the  velocity  solution,  there  still  is  a  shear  present  at  the 
bottom  of  the  mixed  layer  due  to  the  inertial  component  of  the  velocity  solutions. 
Once  an  inertial  wave  is  excited,  it  will  exist  till  it  is  damped  out  by  the  ambi¬ 
ent  ocean.  An  average  damping  time  is  about  10  days  (Pollard,  1970).  At  this 
point,  one  might  tend  to  think  that  model -simulated  shear  exists  only  at  the  top 
(surface)  and  bottom  of  mixed  layers.  This  is  not  the  case,  as  is  shown  in  the 
next  example. 

On  day  number  18.58,  the  wind  stress  began  increasing  after  a  period 
when  it  was  constant,  and  then  began  decreasing.  The  largest  shear  is  now  located 
around  the  penetration  depth  of  the  Ekman  velocity,  which  is  about  10  meters  (fig- 
rue  20).  The  next  largest  shear  is  past  30  meters  and  is  caused  by  a  previously- 
excited  inertial  wave. 

The  qualitative  picture  is  one  of  an  Ekman  velocity  being  matched  to 
a  previously-excited  inertial  wave  (figure  19).  The  picture  can  be  pursued  fur¬ 
ther  by  analyzing  the  excitations  of  inertial  waves  in  time  during  MILE.  Figures 
21  through  22  show  the  MILE  velocity  hodographs  at  depths  of  1  and  30  meters.  The 
graphics  do  not  have  enough  resolution  in  time,  and,  therefore,  lines  that  should 
be  arcs  of  a  circle  appear  as  straight  lines.  In  figure  21,  the  time  span  from 
2  to  7  days  is  covered  and  a  range  of  inertial  amplitudes,  indicating  that  as  wind 
stress  magnitude  varies  the  amplitudes  of  the  excited  inertial  waves  vary.  In 
figure  22,  the  time  span  from  32  to  35  days  at  a  depth  of  30  meters  is  covered, 
and  one  can  see  that  a  constant  amplitude  inertial  wave  was  present.  This  iner¬ 
tial  wave  was  excited  by  a  large  storm  from  day  22  to  day  25.  Since,  at  that 

time,  the  temperature  mixed  layer  was  already  pre-mi xed  to  about  30  meters,  the 

o 

large  wind  stress  that  peaked  at  5  dynes/cm  was  able  to  excite  a  large  ampli¬ 
tude  inertial  wave  at  a  depth  of  30  meters.  Similarly,  on  day  number  18.58,  there 
was  a  large  inertial  wave  present  that  had  been  excited  during  the  wind  storms 
from  day  9  to  day  14.  This  inertial  wave  penetrated  to  a  depth  of  almost  30 
meters  (figure  19). 

IV.  CONCLUSION 

The  shear  was  forecasted  for  the  MILE  experiment  with  a  one-dimensional 
higher-order  closure  mixed  layer  model  by  Warn-Varnas  and  Piacsek  (1979).  The 
statistical  distributions  of  shear  and  stratification  were  analyzed  in  various 
ways.  The  probability  distribution  of  flux  Richardson  numbers  showed  a  shift 
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toward  higher  Richardson  numbers  as  the  depth  of  penetration  into  the  ocean  in¬ 
creased.  This  situation  arises  because  the  depth  of  penetration  of  shear-produced 
turbulence  into  the  ocean  depends  upon  the  strength  of  the  wind  stress.  The 
larger  values  of  wind  stress  are  less  probable  and,  therefore,  the  probable  values 
of  the  Richardson  number  increase  with  depth,  indicating  that,  on  the  average, 
there  is  less  turbulence  with  increased  depth. 

The  distribution  of  shear  and  stratification  during  MILE  were  obtained  by 

computing  scatter  diagrams  at  various  depths.  These  scatter  diagrams  showed  that, 

most  of  the  time,  the  largest  shears  are  located  around  the  bottom  of  the  mixed 

layer  and  near  the  surface.  Within  the  mixed  layer,  the  shears  are  usually 

smaller,  indicating  a  high  homogeneity  of  velocity  profiles.  For  transient  cases, 

there  are  exceptions,  as  will  be  pointed  out  later.  The  maximum  shears  were  dis- 

-2  -2 

tributed  in  the  range  of  2  x  10  to  7  x  10  (cm/sec).  These  values  are  in  the 
range  of  the  experimental  results  of  R.  E.  Davis,  et  al.  (1980). 

Steady  and  unsteady  cases  of  high  and  low  wind  stress  were  analyzed.  The 
analysis  was  performed  by  decomposing  the  velocity  into  an  Ekman  and  inertial 
part.  The  Ekman  or  shear  part  was  dependent  upon  the  wind  stress  being  present 
at  the  surface,  while  the  inertial  part  was  influenced  by  present  and  past  wind 
stresses.  Once  an  inertial  wave  was  excited,  it  remained  for  about  10  days  till 
it  was  damped  out.  The  hodographs  showed  that  as  the  wind  stress  amplitudes 
changed,  the  amplitudes  of  excited  inertial  waves  also  changed.  The  Ekman  part  of 
the  velocity  diminished  to  zero  as  the  wind  stress  decreased  to  zero. 

A  transient  case  was  found  in  which  the  maximum  shear  does  not  occur  at  the 
bottom  of  the  temperature  mixed-layer  depth.  This  can  happen  when  there  is  a 
building  Ekman  part  of  the  velocity  solution,  generated  by  a  rising  wind  stress, 
being  superposed  on  a  previously-excited  inertial  wave. 

V.  ACKNOWLEDGMENT 

This  work  was  sponsored  by  NORDA  Code  500. 

VI.  REFERENCES 

Davis,  R.  E.,  deSzoeke,  R.,  Halpern,  D.,  and  Niiler,  P.,  "Variability  and 
Dynamics  of  the  Upper  Ocean  During  MILE,"  submitted  to  Journal  of 
Physical  Oreanooraphy  1980. 

Mellor,  G.  L.,  "Analytic  prediction  of  the  properties  of  stratified  plane¬ 
tary  surface  layers,"  J.  Atmos.  Sci . ,  30,  1061  (  1973). 


8 


Miyake,  M.,  "Oceanographic  observations  at  Ocean  Station  P,"  Pacific  Marine 
Science  Report  78-1,  Institute  of  Ocean  Sciences,  Sidney,  B.  C.  (1978). 

Nomitsu,  T.,  "A  theory  of  the  rising  stage  of  drift  current  in  the  ocean," 
Mem.  Coll.  Sci . ,  Kyoto,  16,  No.  2,  (1933). 

Pollard,  R.  T.,  and  Millard,  R.  C.,  "Comparison  between  observed  and  simu¬ 
lated  wind-generated  inertial  oscillations,"  Deep-Sea  Research  17, 
813-821,  (1970). 

Warn-Varnas,  A.  C.,  Dawson,  G.  D.,  and  Martin,  P.  T.,  "Forecast  and  Studies 
in  the  Oceanic  Mixed  Layer  During  the  MILE  Experiment,"  Geophys. 
Astrophs.  Fluid  Dynamics.  To  appear  in  February  issue  1981. 

Warn-Varnas,  A.  C.  and  Piacsek,  S.  A.,  "An  investigation  of  the  importance 
of  third-order  correlations  and  choice  of  length  scale  in  mixed  layer 
modelling,"  Geophys.  Astrophys.  Fluid  Dynamics,  13,  225  (1979). 


9 


W  IND  STRKSS 


0  0.0  10.0  III 

0  12.0  14 

.0  140  15. 

0  16.0  |7.0 

I'OTAL  H.UXKS 

.0  0.0  10.0  1 1.0  121 
A 

)  14.0  14.0  15 

.0  16.0  17.0 

WIND  STRKSS 

2 

>  20 


184 

l« 

I4( 

>  1200.0 
£1000.0 

s  600.0 

S'  4( 


17.0  IH.0  10.0  20.(1  21.0  22.0  2.1.0  24.0  25.0  2 6.0  27.0  2K.0  20.0  40.0  41  jO  42.0 

TOTAL  Kl.llXt:s 


17.0  1X  0  10.0  20.0  210  22  0  2  VO  2  (.0  25.0  26.0  27.0  2S.0  2‘>.0  40.0  INI  42.0 

B 

Figure  1.  The  atmospheric  forcing  fluxes  of  MILE. 
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Figure  3.  Distribution  of  flux  Richardson  numbers  at  10  meters. 
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Figure  14.  Total  velocities  versus  depth  on  day  7.92. 
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Figure  22.  Velocity  hodographs  for  MILE  from  day  32  to  36  at  a  depth  of  30  meters. 
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APPENDIX 

THE  FINITE-DIFFERENCE  SCHEMES 

For  the  WP  model,  the  finite-difference  equations  are  solved  implicitly 
for  un+1,  vn+1,  Tn+*.  This  involves,  as  explained  in  the  WP  paper,  solving  alge¬ 
braically  the  finite-difference  equations  for  the  second-order  correlations.  The 
resultant  finite-difference  equations  for  the  mean  fields  are,  then,  in  the  desir¬ 
able  diffusive-type  form  plus  a  forcing  term.  The  diffusive  form  is  desirable 
because  it  enables  one  to  solve  implicitly  for  the  mean  fields.  As  an  illustra¬ 
tion,  the  resultant  equation  for  temperature  is  as  follows: 
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The  treatment  of  the  Coriolis  terms  are  changed  from  the  one  stated  in  the 
WP  paper.  The  method  of  time  splitting  is  now  used  to  separate  them  from  the  rest 
of  the  solution.  This  treatment  is  illustrated  in  equation  (2.1).  The  finite- 
difference  analogy  is 
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where  is  the  eddy  coefficient  term,  is  the  source  term,  and  u*  is  the  inte  . 

mediate  time  value.  This  method  was  introduced  so  as  not  to  overdamp  the  inertial 
oscillations. 

The  grid  is  stretched.  The  stretching  of  the  grid  and  the  resultant  accu-  i. 

racy  of  the  overall  finite-difference  scheme  is  described  in  the  WP  paper. 
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